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ABSTRACT 

Aims. We present a study of the Be/X-ray binary system V 0332+53 with the main goal of characterizing its behavior 
mainly during the intermediate-luminosity X-ray event on 2008. In addition, we aim to contribute to the understanding 
of the global behavior of the donor companion by including optical data from our dedicated campaign starting on 2006. 
Methods. V 0332+53 was observed by RXTE and Swift during the decay of the intermediate-luminosity X-ray outburst 
of 2008, as well as with Suzaku before the rising of the third normal outburst of the 2010 series. In addition, we present 
recent data from the Spanish ground-based astronomical observatories of El Teide (Tenerife), Roque de los Muchachos 
(La Palma), and Sierra Nevada (Granada), and since 2006 from the Turkish TUBITAK National Observatory (Antalya). 
We have performed temporal analyses to investigate the transient behavior of this system during several outbursts. 
Results. Our optical study revealed that continuous mass ejection episodes from the Be star have been taking place 
since 2006 and another one is currently ongoing. The broad-band 1-60 keV X-ray spectrum of the neutron star during 
the decay of the 2008 outburst was well fitted with standard phenomenological models, enhanced by an absorption 
feature of unknown origin at about lOkeV and a narrow iron K-alpha fluorescence line at 6.4 keV. For the first time in 
V 0332+53 we tentatively see an increase of the cyclotron line energy with increasing flux (although further and more 
sensitive observations are needed to confirm this). Regarding the fast aperiodic variability, we detect a Quasi-Periodic 
Oscillation (QPO) at 227+9 mHz only during the lowest luminosities. The latter fact might indicate that the inner 
regions surrounding the magnetosphere are more visible during the lowest flux states. 
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1. Introduction 

Accreting X-ray pulsars are binary systems composed of 
a donor star and an accreting neutron star (NS). In High 
Mass X-Ray Binary (HMXB) systems the optical compan¬ 
ion could be either a massive early-type supergiant (su¬ 
pergiant systems) or an 0,B main sequence or giant star 
(BeX binaries; BeXB). Among the most remarkable signa¬ 
tures found in BeXBs are the detection of IR excess and 
emission-line features in their optical spectra produced in 
a disc-like outflow around the Be star. Historically, their 
outbursts have been divided into two classes. Type I (or 
normal) outbursts normally peak at or close to periastron 
passage of the NS {Lx < 10^”^erg s“^). Type B (or giant) 
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outbursts reach luminosities of the order of the E ddington 
luminosity (Lx ^10^® erg s~^: iFrank et al.l [200211 . with no 
preferred orbital phase. We will refer as “intermediate” to 
any X-ray outburst that does not comply with the afore¬ 
mentioned properties, thus has (Lx erg s“^). 

1.1. V 0332+53 

The recurrent hard X-ray transient V 0332+53 was de¬ 
tected with the Vela 5B observatory in 1973 during its giant 
outburst, reaching a peak intensity of »1.6 Crab in the 3- 
12 keV energy band (iTerrell fc PriedhorskyLfTo^ . The sys¬ 
tem had passed a ten-year X-ray quiescent phase when 4.4 s 
pulsa tions were detected with Tenma and EXOSAT satel¬ 
lites (IStella et al.l . I1985D . These X-ray activities, a series 
of Type I outbursts, lasted about three months separated 
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by the orbital period (34.25 d) of the system. During these 
outbursts, rapid random fluctuations in the X-ray emis¬ 
sion in a ddition to the puls e-profile variations were also 
reported (lUnEfer et al.l . Il992ll . The system underwent an¬ 
other outburst, classified as Type II, in 1989 which led to 
the discovery of a cyclotron line scattering feature (CRSF) 
at 28.5 keV and two QPOs centred at 0.05 1 Hz and 0.22 Hz 
( Makishima et al.l . 119901 : iTakeshima et al.l . 11994 iQu et all . 
[ 200 ^ 


The optic al _ companion of the system, BQ Cam 

(jArgvle et al.l . Il983^ is a n 08-9 type main , seque nce star 
at a distance of ~7kpc (|Negueruela et al.l . Il999ll . It has 
been widely observed both in optical and IR wavelengths 
since its identification. The optical spectrum is charac¬ 
terised by the highly variable emissions of Hq,, H,g and 
H.y in addition to the He I lines. The photometric data 
shows an IR excess (l^e rnacca et al. . Il984 Coe et all . 
ll987t[Honevcutt fc Schlege]Lll985l : llJnger et al.Lll992tl . The 
brightening in optical/IR light-curves is usually accompa¬ 
nied with the X-ray outburst pha ses as in the case of gi¬ 
ant 20 04 outburst of the system. iGoranskii fc Barsukoval 
(|2004ll predicted this outburst based on the optical bright¬ 
ening of V 0332+53 in optical/IR band. During the out¬ 
burst, three additional CRSFs at 27, 51 and 74keV were 
detected in Rossi X-ray Timin g Explorer (RXTE) ob¬ 
servations (jCoburn et al.l . l2005ll and confirmed by th e 
subsequent INTEGRAL data (|Pottschmidt et al.l . l2005ll . 
iTsvgankov et al.l (j2006ll showed that the energy of these 
features are linearly anti-correlated with its luminosity. The 
following X-ray activities of the system were in 2008, 2009 
and 2010 w i th relatively weaker pe ak fluxes (iKrimm et al.l . 
l2008l 120091 : iNakaiima et al.L l2010ll . The system was in X- 
ray quiescence until 18 June 2015, when the Be companion 
probably reached its ma ximum optical brightness aft er re¬ 
newing activity in 2012 (ICamero-Arranz et al.l . 1^151) . 

Here, we present a multi-wavelength study of 
V 0332+53 mostly during the recent events initiated in 
2008. For this purpose we used archived Swift-IEBU and 
RXTE pointed observations carried out in 2008, as well as 
one Suzaku observation from 2010 and survey data from 
different space-borne telescopes and covering intermedi¬ 
ate and low-luminosity events. In addition, we used opti- 
cal/IR data from our dedicated campaign involving several 
ground-based astronomical observatories. This provides us 
a unique opportunity to undertake studies of X-ray out¬ 
bursts of intermediate and low X-ray luminosities. The lat¬ 
ter is a regime very scarcely explored in V 0332+53, one of 
the BeX with strongest magnetic field already known. 


2007 2008 2009 2010 2011 2012 2013 2014 2015 



Fig. 1. Comparison of the ^wi/it/BAT light curve (15- 
50keV) with the optical magnitudes (ROTSE, OMC, 
IAC80) of V 0332+53 since 2006. The BAT, ROTSE and 
OMC data points have been rebinned for clarity. The vi¬ 
olet arrows in the X-ray panel represent the times of the 
periastron passages of the NS. 


aperture photometry (ICamero et al.L I20l4 for more details 
on reduction). 

The main part of the long-term optical CCD ob¬ 
servations of V 0332+53 include data from the 0.45 m 
ROTSEHIc0 telescope from the TUBItAK National 
Observatory (TUG; Antalya, Turkey) between February 
2006 until December 2012. This telescope operates with¬ 
out filters and is equipped with a 2048 x 2048 pixel CCD. 
Dark and flat-field corrections of all images were done au¬ 
tomatically by a pipeline. Instrumental magnitudes of all 
the corre cted images were obtain ed using the SExtractor 
Package (iBertin fc Arnout^ . 1 199611 . Calibrated ROTSEHId 
magnitudes were acquired by comparing all the stars in 
each frames with USNO-A2.0 catalog i?-band magnitudes. 

In addition, we have us ed data from the Op tical 
Monitoring Camera (OMC; iMas-Hes se et al. l l2003ll o n 


board the INTEGRAL satellite (jWinkler et al. . 200311 . 


They were obtained from the Data Archive Unit at CAB 
Astrobiology centr^^ (INTA-CSIC), and pre-processed by 
the INTEGRAL Science Data Centre (ISDCj^. 


2. Observations and Data Reduction 

2.1. Optical/IR photometric observations 

The optical counterpart has been observed in the optical 
band with the 80-cm IAC80 telescope at the Observatorio 
del Teide (OT; Tenerife, Spain) starting on August 2014 
until February 2015 (see on-line material Tab.[T]and Fig.[T]). 
We obtained the optical photometric CCD images using B 
and V filters with integration time of 60 s. In addition, on 
20 December 2014 a single infrared J, H and Kg observation 
was performed by the CAIN camera from the 1.5 m TCS 
telescope at the OT, with integration time of 150, 90 and 
90s, respectively. The reduction of the data was done by 
using the pipelines of both telescopes based on the standard 


2.2. Optical spectroscopic observations 

Optical spectroscopic observations of the donor star were 
performed since 26 September 2006 until 22 November 
2014 (see on-line material TabHJ with three different tele¬ 
scopes: the Russian-Turkish 1.5 m telescope (RTT150) at 
the TUG, the 2.56 m Nordic Optical Telescope (NOT) lo¬ 
cated at the Observatorio del Roque de los Muchachos (La 


^ The Robotic Optical Transient Search Experiment, ROTSE, 
is a collaboration of Lawrence Livermore National Lab, 
Los Alamos National Lab, and the University of Michigan 
(http: //www. ROTSE.net) 

^ http://sdc.cab.inta-csic.es/omc/index.jsp 
® http://www.isdc.unige.ch 
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Table 1. Log of the X-ray observations^. 


Obs. 

Date ^ 

MJD 

Telescope 

Exp. 

Pulse 





time 

period 





(s) 

(s) 

1 

04/11/08 

54774.636 

RXTE^ 

1600 

4.3742(1) 

2 

07/11/08 

54777.999 

RXTE 

752 

- 

3 

09/11/08 

54779.731 

RXTE 

496 

- 

4 

10/11/08 

54780.594 

RXTE 

864 

- 

5 

12/11/08 

54782.394 

RXTE 

1456 

4.3740(1) 

6 

14/11/08 

54784.511 

RXTE 

1552 

4.3746(6) 

7 

16/11/08 

54786.543 

RXTE 

1488 

4.37450(2) 

8 

19/11/08 

54789.826 

RXTE 

3 840 

4.37540(2) 

9 

20/11/08 

54790.775 

RXTE 

2 848 

4.3763(1) 

10 

21/11/08 

54791.743 

RXTE 

1152 

4.3760(1) 

11 

22/11/08 

54792.270 

RXTE 

3 716 

4.37610(4) 

12 

23/11/08 

54793.273 

RXTE 

1597 

4.37620(4) 

13 

23/11/08 

54793.387 

RXTE 

3 235 

4.37620(5) 

(5) 

12/11/08 

54782.384 

SwifC 

1706 

- 

(6) 

14/11/08 

54784.520 

Swift 

2108 

- 

(7) 

16/11/08 

54786.533 

Swift 

2 307 

- 

(8) 

19/11/08 

54789.338 

Swift 

1613 

- 

(9) 

20/11/08 

54790.811 

Swift 

1597 

- 

(10) 

21/11/08 

54791.487 

Swift 

1489 

- 

14 

16/02/10 

55243.264 

Suzaku ® 

15 946 

4.37(1) 


^ Standard mode for RXTE and Suzaku and Windowing 
mode for Swift-XHT. 

^ Observation date in the format (dd/mm/(20)yy). 

® RXTE/PCA exposure time. 

Swift/XRT exposure time. 

® Suzaku/XlSl exposure time. 

Palma, Spain), and the 1.5 m Telescope at the Observatorio 
de Sierra Nevada (OSN-CSIC) in Granada (Spain). 

The spectroscopic data from RTT150 were obtained 
with the TUBiTAK Faint Object Spectrometer and 
Camera (TFOSC). The reduction of the spectra was done 
using the Long-Slit package of MIDA^. The low-resolution 
OSN spectra were acquired using Albireo spectrograph 
whereas the NOT spectra were obtained with the Andalucia 
Faint Object Spectrograph and Camera (ALFOSC). The 
reduction of this data set was performed usin g standard 
procedures within IRAP0 (jCamero et al.l 1^014 for further 
details on the instrumentation and data reduction). 

All spectroscopic data were normalized with a spline 
fit to continuum and corrected to the barycenter after the 
wavelength calibration. The full width at half maximum 
(FWHM) and equivalent width (EW) measurements of Hq, 
lines were acquired by fitting Gaussian functions to the 
emission profiles using the ALICE subroutine of MIDAS. 

2.3. X-ray Observations 

The XRT (on-board the Swift satellite: ICehrels et al.ll2003) 
observed V 0332+53 in Windowing Timing (WT) mode 
(1.7ms time resolution, ID Image) six times from 12 to 
21 November 2008. Tab. [1] shows a summary of the log of 

http://www.eso.org/projects/esomidas 
® IRAF is distributed by the National Optical Astronomy 
Observatory, optical images which is operated by the 
Association of Universities for Research in Astronomy (AURA) 
under cooperative agreement with the National Science 
Foundation. 


these observations. To reduce the data and extract final 
products (spectra and light curves) we followed the stan¬ 
dard recipe as detailed in the Swift/XBT data reduction 
guid43. The relevant response matrix to use is given by the 
HEASARC’s calibration database (CALDEQ). All the XRT 
lightcurves were barycentred corrected. We also used for 
comparison data products supplied by the UK Swift Science 
Data Centre at the University of Leicester (lEvans et al.L 

IMl)- 

_ RXTE Proportional Counter Array fPCA: [jahoda et al.l 

1199611 standard 1 and 2 data only from PCU-2 were selected 
for this work. To extract PCA products we followed the 
standard procedure showed in its cookboolfl. A systematic 
error of 0.5% was added in quadrature to the PCA stan¬ 
dard 2 spectra, as recommended by the instrument team. 
Products from the High Ener gy X-ray Timing Experiment 
(HEXTE; iGruber et al.lll99^ were also extracted using the 
corresponding cookbook recip^B- Only standard mode data 
from cluster B were used. All the RXTE lightcurves were 
barycentred corrected. 

In ad dition, we analysed o ne observation performed by 
Suzaku (|Mitsuda et al.l . l2007ll in 16 February 2010. We 
select ed only data from th e X-ray Imaging Spectrometer 
(XIS; iKovaina et al.l l2007ll CCD and the PIN diode de- 
tector (10-70 keV) of the he Hard X-ray Detector (HXD; 
iTakahashi et al.l l2007ll . For this observation only XISO, 
XISl and XIS3 detectors were available. Following the 
Suzaku ABC auidj^ we reprocessed all the data and ex¬ 
tracted the final products. The clocking mode option was 
normal burst (1.0) with the 1/4 window selected (Is time 
binning), and a data editing mode of 3x3. No pile-up cor¬ 
rection was needed. The events were barycentred c orrected . 

W e used the HEAsoft software version 6.15.1 (|Arnaudl . 
Il996h for the data analysis of all the instruments. Since 
2008, the Gamma-ray Burst Monitor (CBM) on board the 
Fermi satellite, has been monitoring V 0332+53. In this 
work we used tin iing products provided by the CBM Pulsar 
TeanJ^ (see e.g. iFinger et al.l . l200^ ICamero-Arranz et al.L 
I 2 OIOL for a detailed description of the timing technique). We 
also used quick-look X-ray results provided by the RXTE 
All Sky Monitor teairP^. and since 2008 we used the quick- 
look X-ray Swift/BAT tra nsient monitor resu lts provided 
by the Swift/BAT teanO (jKrimm et al.L[20T^ . See Fig.[2l 

3. Results 

3.1. Optical photometry 

The optical light-curves of BQ Cam compared to the 
Swift/BAT data (15-50 keV) from June 2006 to June 2015 
are shown in Fig. [TJ In 2008, after being w3.5yrs in qui¬ 
escence in X-rays the system underwent a new active pe¬ 
riod. However the optical companion entered the brighten¬ 
ing state nearly 1.5 yr before the NS. When the enhance¬ 
ment in the optical brightness reached a value of RiO.25 mag. 

® http://swift.gsfc.nasa.gov/analysis/xrt_swguide_vl_2.pdf 
’’ http://heasarc.gsfc.nasa.gov/FTP/caldb 
® https://heasarc.gsfc.nasa.gov/docs/xte/recipes/ 
pca_spectra.html 

® https://heasarc.gsfc.nasa.gov/docs/xte/recipes/hexte.html 
http://heasarc.nasa.gov/docs/suzaku/analysis/abc 
http://gammaray.nsstc.nasa.gov/gbm/science/pulsars 
http://xte.mit.edu/ASMJc.html 
http://swift.gsfc.nasa.gov/results/transients 
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Fig. 2. Top. Long-term pulse period history of V 0332+53 since its discovery in 1973. Middle. Long-term optical 
lightcurve of V 0332+53. All but the IAC80 light curve have been rebinned for clarity. Bottom. Type of X-ray outburst 
underwent by V 0332+53 (i.e. Type II, Intermediate, and Type I). The intensities are the same for outbursts of the same 
type, not having been corrected for the different energy bands. We only intend to be illustrative of the times and types 
of events. The long-term evolution o f the EW (red fil l ed circles) is also ov e rplotted. The measurem e nts have been 


extrac t ed from live fc Kodaird ( 198^. Kodaira et al. (ll985ll . [Corbet et H] (ll986D . lNegueruela et al.l (ll999D . lMasetti et al.l 
(j2005ll , iKiziloglu et al.l ( 200811 . Goranskii et al.l ( 201011 and the present work. 


The X-ray outburst was triggered, around 17 October 2008 
(MJD 54756) roughly six days before the periastron passage 
of the NS (vertical red-dashed lines in Fig. [IJ. This large 
flare lasted about 40 d and reached a flux of R:!214mCrab 
within 3 weeks. It should be pointed out that the optical 
magnitude of V 0332+53 shows a sudden decrease when the 
X-ray outburst declines. Fading of the optical companion 
continued until the end of January 2009 (MJD «54850) and 
soon after it turned back to its brightening state. The X-ray 
flux, on the other hand, stayed at its quiescent level, until 
the Be star reached a peak value of «:!l3.8mag (ROTSE 
magnitude). A new X-ray phase of the system started in 
November 2009 again a few days after the periastron pas¬ 
sage (MJD r::! 55140) of the NS. That activity included five 
small but periodic outbursts separated by the orbital period 
(34.25 d) of the system classified as Type 1. 

The second outburst of these series was the most power¬ 
ful one with a peak flux of ~110 mCrab. The optical magni¬ 
tude, on the other hand, started to decrease with declining 
of the second outburst (MJD ss55201). The X-ray activi¬ 
ties finished by the end of the May 2010 (MJD Ri55312) 
while the optical magnitude was still fading. It ceased 
around August 2011 (MJD Ri55780) and lasted ~4 years. 
Our recent data from the IAC80 telescope and the opti¬ 
cal camera OMC on board INTEGRAL, shows an optical 
enhancement of the counterpart BQ Cam, with the bright¬ 
ness varying from the quiescence level of V = 15.93±0.55 
(06 Jan. 2012; MJD 55932. 60) to V = 15.332A0 . 008 ( 1 
Feb. 2015; MJD 57054.04; ICamero-Arranz et al.l [20l|). 
This corresponds to an increment of ssO.Omag, simi¬ 


lar to that observed durii ig the brightenin g episodes oc- 
curred e.g. in 1983 (se e GoranskiT l200lD and in 2004 


(|Goransk ii fc Ba rsukova . 2004 ). but high er than the ones 


m 2008 (|Kaur et al.l . 120081) and 2009 (iGoranskii et al.L 
120101) . This is confirmed by our only IR measurement 
from December 2014 (MJD 57011.869), that showed an 
unusually bright Be companion. The magnitudes ob¬ 
tained were J=11.228±0.003mag, 11=10.651±0.003mag, 
and Ks=10.285±0.00 4mag, which are comp arable to those 
from December 1983 (|Williams et al.l . Il983l) . We note that 
BQ Cam did not reach its maximum level (V = 15.2 mag) 
in our observations from February 2015. The maximum has 
been probably reached recently, since the BAT monitor is 
currently detecting X-ray activity from V 0332+53 with a 
daily average flux of the order of 310.4 mCrab on 25 June 
2015 (MJD 57198.0). 


3.2. H„ line 

The long-term Hq line monitoring of BQ Cam from 
September 2006 to November 2014 has been studied (see 
Fig. 131). In contras t to the previous works statin g the Hq 
profile variations (iNegueruela et al.l . I1998L Il999l) . we did 
not see such variability patterns during the observation pe¬ 
riod. Instead, all the line profiles are nearly-symmetric in a 
single-peaked form despite the lack of the observations for 
the period 2010-2011. 

In spite of having the same single-peaked emission pro¬ 
file on a long-term basis, the Hq line strength does not show 
a constant trend during the spectroscopic observations (see 
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Fig. 3. Hq, profiles of V 0332+53 selected between 2006 and 2014. The single-peaked emission of the line does not change 
on a long-term basis. 


on-line material Tab. [a. In Fig. 0] the evolution of EW and 
FWHM measurements of the emission line are shown. The 
Hrst two measurements of the line coincide with the declin¬ 
ing phase of the optical outburst seen in the time interval 
MJD 53576-54180. The EW values decreased as the star 
faded away which can be attributed to the weakening of the 
decretion disc. Besides, this decreasing pattern continues af¬ 
ter the optical outburst triggered (MJD 54324). Although 
our data coverage is not uniform, we see a reverse relation 
between the optical magnitude and the width of the emis¬ 
sion line until 4 March 2009 (MJD 54894) when the lowest 
value of the EW of Hq, is reached (see Fig. [21). After about 
seven months, the EW reached a value of -5.21 A. Although 
this sudden increase in the line width well fits with the ris¬ 
ing of the optical magnitude, we cannot make any further 
explanation regarding the behavior of the emission line due 
to the lack of data for 2010-2011. In 2012-2014, we see that 
the strength of the Hq, line kept its varying pattern without 
a noticeable change in the optical magnitude indicating the 
upcoming mass ejection event. 

In contrast to the variations in EW values, the FWHM 
values stayed approximately constant until the end of 2012. 
At that point, the FWHM values started to follow nearly 
the same pattern of EW. Indeed, showing a similar behav¬ 
ior is not unusual for EW and FWHM, since the width of 
the emission profile is directly related to the rotational ve¬ 


locity of the Be star. iHanuschil] (I1989H gives the relation 
between the Hq widths and the projected rotational veloc¬ 
ity, usin(f). The projected rotational velocity of V 0332+53 
is found to be t;sin(z) = 166.68kms“^, using the average 
values of-5.04A (i.e. 230kms“^) and 7.3 A (333.56kms“^) 
for the EW and FWHM, respectively (see on-line material 
Tab. |2|). For the true rotational velocity of V 0332+53, we 
have assumed th e inclination a ngle o f the system of i< 18.9° 
as suggested bv IZhang et al.l (|2005ll . Therefore, the lower 
limit to the true rotational velocity would be 514kms“^. 
Using the typical mass and radius values for a late 0-type 
star (M>20Mq and R>15R0), the break-up velocity is 
determined as 600kms“^ that is close to the rotational ve¬ 
locity. 

3.3. X-ray activity 

3.3.1. Periodic timing analysis 

To search for periodic pulsations from the compact ob¬ 
ject we have used the tool efsearch from the Xrono0 
package. This tool searches for periodicities in a time se¬ 
ries by folding the data over a range of periods and by 


http://heasarc.gsfc.nasa.gov/docs/xanadu/xronos/ 
xronos.html 
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2007 2008 2009 2010 2011 2012 2013 2014 2015 



Fig. 4. Evolution since 2006 of the EW and FWHM mea¬ 
surements of the Hq, emission line. The red-dashed verti¬ 
cal line denotes the beginning and the end of the optical 
outburst seen in ROTSEIIId light curve (see also Fig. [T|) 
whereas the black-dashed line shows the time of the turning 
point from decreasing to increasing trend just after (RilOd) 
the strong X-ray outburst around MJDRi54850. 

searching for a maximum chi-square as a function of pe¬ 
riod. Standard 1 PCA, XRT and Suzaku-XlS barycentred 
lightcurves for V 0332+53 were selected for this task (see 
Tab. [T]). The resulted pulse periods can be seen in Fig. [2] 
and Tab. 1. The error on the measured spin periods is calcu¬ 
lated using the chi-square versus spin period plot provided 
by efsearch. The points near the peak of this plot are 
fitted with a gaussian. The error on the gaussian centre es¬ 
timate is taken as the error on the spin period. From Tab.[T] 
we can see that during the 2008 outburst decay the neu¬ 
tron star seems to spin slower, as expected. Additionally, 
we note that in Fig. [5] we do not include the pulse period 
obtained with Suzaku-^IS (Ppuise =4.37(1) s) because of 
its large uncertainty. Furthermore, we do not include any 
pulse period from the XRT data because they were not 
significantly detected. 

3.3.2. Aperiodic timing analysis 

We performed an analysis of the fast time variabil¬ 
ity of V 0332+53 in the 2-60 keV energy range of 
Obs. 1, 2, 6, 8 and 9. The net PCA count rate is 
Ri700,600,300,300,160ctss“^ for Obs. 1-9, respectively. 
The rest of the observations either are too short and/or 
the source is too faint to perform timing analysis. The time 
resolution of the PCA in the mode used is 0.125 s. 

We used the GHATS package, developed under the 
IDL environment at INAF-OAEPH. to produce the Power 
Density Spectra (PDS) from 512 points in each light curve. 
The PDS were then averaged together for each observa¬ 
tion. The PCA light curve was binned at its time resolution 
(0.125s). This yields a Nyquist frequency of = 4Hz. The 

http://www.brera.inaf.it/utenti/belloni/ 
GHATS-Package/Home.html 


PDS were normalized according to iLeahv et al.l (Il983ll . All 
of the PDS show low-frequency band-limited noise (that 
may turn into peaked noise) and QPO noise. 

PDS fitting was carried out with the standard XSPEC 
fitting package by using a unit response. Fitting the (2- 
60keV) PDS with a model constituted by a zero centred 
Lorentzian for the flat-topped noise (hereafter called Li) 
plus a constant for the Poissonian noise (i.e. well repro¬ 
duced by a powerlaw model component with its photon 
index frozen to zero in the fits) results in mostly accept¬ 
able chi-square values only for the PDS of Obs. 1, (see 
Tab. HI). In the case of Obs. 8, 9 the fit statistics with 
this model is worse and leaves some visible residuals, i.e. 
X^/v = 290/252, 233/252, respectively. The fit to these ob¬ 
servations improved by adding a further Lorentzian com¬ 
ponent (L 2 ) centred at 1 / = 30 — 80 mHz. This compo¬ 
nent changed the (F-test) fits statistics by Ax^Ri20 — 30 
for 1 / = 249 d.o.f., thus a <3 (t improvement. For Obs. 2, 6 
we took into account the latter feature only, since it is broad 
and therefore constitute a main part of the measure of the 
rms, i.e. by (10 — 20)%. This resulted in a ht statistics of 
^ jv = 253/252,255/252 for Obs. 2, 6, respectively. 

The ht of the (2-60 keV) PDS with a model consti¬ 
tuted by a zero (Li) and a ^ = 30 — 80mHz (L 2 ) centred 
Lorentzians (L 2 only for Obs. 2 and 6) for the band-limited 
noise plus a constant for the Poissonian noise results in a 
poor description of the data for Obs. 6, 8, 9. The Ht statis¬ 
tics with this model is of y^/j^ = 255/252,263/249, 209/249 
for Obs. 6-9, respectively. In these observations there are 
positive residuals centred in the range of (220 — 240) mHz 
that we Htted by adding a further Lorentzian compo¬ 
nent (Lqpo) centred at those frequencies. This compo¬ 
nent changed the Hts statistics by Ax^Ril2,19,13 for v = 
249,246,246 d.o.f., thus a (2.9,3.9,3.0)cr improvement for 
Obs. 6-9, respectively. We took into account this feature, 
since it is broad and therefore affects substantially the mea¬ 
sure of the rms, i.e. by Ri(2 — 4)% for all the observations. 
These positive features have a peaked form, thus we will call 
theni QPO s , consistently with previous Hndings (|Ou et al.L 
120051 : iReid . 120081 : iReig fc Nesooli 1201311 . The detection of 
this peak is only signihcant for Obs. 8, which can be ex¬ 
plained since this is the observation with the longest expo¬ 
sure time of the sample. 

We calculated the fractional rms from the best Ht 
model (integrated in the (8x10“^ — 4) Hz band). This was 
found to be (15 — 20)% in the 2-60keV energy range. 
The results obtained from the global fits (i.e. using the 
lorentz+lorentz+lorentz+powerlaw model in XSPEC) 
are reported in Tab. HI We plot in Fig. [5] the broad-band 
PDS with the best-fit model. We notice that our results 
are broadly in agreem ent with those previously obtained 
()Reig fc NesDolil . l2013ll . 

3.3.3. Spectral analysis 

We fitted the brightest background-subtrac ted spectra wit h 
standard spectral models in XSPEC 12.8.2 (|Arnaudl . 119961) . 
All errors quoted in this work are 68% (Icr) confidence. 
The spectral fits were limited to the 1-10, 4.5-28, 26-60 keV 
range for the XRT, PCA and the HEXTE, respectively, 
where the calibration of the instruments is the best. The 
spectra were grouped in order to have at least 100, 50, 
50 counts for the XRT, PCA and HEXTE, respectively for 
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Table 2. Results from the RXTE timing analysis^. 



Obs. 1 

Obs. 2 

Obs. 6 

Obs. 8 

Obs. 9 

I'l 

7.5+7 

- 

- 

10+8 

12+10 

FWHM 

174+14 

- 

- 

29+14 

22+20 

Norm. 

21.7+1.2 

- 

- 

1.7+0.8 

1.1+0.6 

rms (%) 

16.4+0.4 

- 

- 

6.2+0.3 

8.3+0.6 


- 

48+13 

36+11 

82+19 

77+25 

FWHM 

- 

147+15 

166+12 

137+25 

160+24 

Norm. 

- 

17.2+1.5 

9.1+0.5 

1.3+0.4 

1.3+0.4 

rms (%) 

- 

16.4+0.7 

16.5+0.5 

6.5+0.3 

8.4+0.6 

t'QPO 

- 

- 

231+5 

226+8 

233+4 

FWHM 

- 

- 

16+11 

33+18 

17+16 

Norm. 

- 

- 

0.46+0.15 

0.16+0.06 

0.12+0.04 

rms (%) 

- 

- 

3.8+2.1 

2.3+1.0 

2.7+1.8 

r 

poisson 

0 

0 

0 

0 

0 

N 

poisson 

1.84+0.04 

1.59+0.05 

1.79+0.03 

1.90+0.03 

1.93+0.03 


0.84(213/252) 

1.0 (253/252) 

0.98 (243/249) 

0.99 (244/246) 

0.80 (196/246) 


^ Model used: lorentz+lorentz+lorentz+powerlaw for five RXTE/VCK observations. 

^ The PDS were created in the (2-60) keV energy band and (8x10“^ — 4) Hz frequency range. 
^ The frequency and width {i/ and FWHM) of the components are shown in units of mHz. 
Errors are 68% confidence errors. 


OBS. 1 (RXTE/PCA) OBS. 2 (RXTE/PCA) OBS. 6 (RXTE/PCA) 




Fig. 5. Power density spectra of the RXTE PCA data in the energy and frequency range (2-60 keV) and 0.008-4 Hz, 
respectively, during observations 1, 2, 6, 8 and 9 (top-left to bottom-right) with the best fit model (solid line) and the 
model components (dotted line). The results of the fitting are shown in Tab. [5] and discussed in the text. 


each background-subtracted spectral channel and to avoid 
oversampling of the intrinsic energy resolution. 

As it is a common procedure in observational studies of 
NS systems we used empirical models to describe the X-ray 
spectra. To compare our results with previous studies we 
used the same (or very similar) spectral components. To fit 
the spectral continuum we used a model composed by pho¬ 
toele ctric absorption (|Bahicinska-Church fc McCammonI 
Il992t phabs in XSPEC) and a power-law with a high-energy 
exponential cutoff (cutoffpl). The absorption component 


was not used in Obs. 1-4 and 11 (we did not fit the data 
from Obs. 12,13 because of the lack of counts), since PCA 
spectra are not affected (i.e. at E>4.5keV). In the case 
of XRT observations, with lower energy coverage (i.e. at 
E>1 keV) the absorption component was included in all our 
fits. This model provided significant wavy residuals in the 
whole spectra of all the observations with non acceptable 
chi-squared values — 60; with i/ being the number 

of degrees of freedom, i.e. d.o.f.). 
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Table 3. Results from the RXTE spectral analysis only^. 



Obs. 1 

Obs. 2 

Obs. 3 

Obs. 4 

Obs. 11 

E(l),.b. 

7.6T2.8 

8.0T2.3 

9.0±0.5 

6.5±0.7 

5±2 

cr(l) 

4.6T1.6 

4.1T0.8 

3.3T0.6 

6±5 

7±3 

S(l) 

15±8 

9±7 

2.2±0.9 

70±40 

llilO 

E(2)gabs 

31.6±0.7 

31.3±0.8 

31.8±0.8 

29.7±1.2 

28.4T1.0 

a(2) 

4.3T0.4 

4.6T0.4 

5.0±0.6 

4.1±0.8 

3.4±0.5 

S(2) 

17±3 

20±4 

23±5 

14±6 

8±3 

■R 

-*-^gauss 

6.3±0.5 

6.40±0.10 

6.9±1.5 

6.3±0.6 

5.7T0.6 

a 

0.18±0.20 

0.22±0.20 

1.6±0.5 

0.5±0.4 

l.liO.6 

N 

0.012T0.005 

0.009±0.004 

- 

0.014±0.010 

0.02±0.01 

r 

4 pOW 

-0.7T0.5 

-0.51T0.24 

-0.93±0.28 

-0.7±0.5 

-0.7T0.3 

EcutofF 

5.7±0.6 

6.6±0.5 

5.7±0.6 

5.0±0.8 

5.68±0.15 

Fx 

10±4 

8±4 

4.7±0.1 

4.5±0.3 

1.3±0.7 

Lx 

5.9±2.4 

4.7±2.4 

2.8±0.3 

2.7±0.3 

0.8±0.4 

CpCA 

1 

1 

1 

1 

1 

Chexte 

0.50±0.05 

0.45±0.07 

1.3±0.3 

1.5±0.3 

0.40T0.07 


0.83 (43/52) 

0.90 (47/52) 

1.1 (56/52) 

1.2 (64/52) 

1.1 (56/52) 


^ Model used: constcuitxgabsxgabs( cutoffpl + gaussiein ) for the 4.5-60 keV RXTE 
spectra. 

^ Centroid, line-width of the lines, normalization of the gaussian line, flux and luminosity 
in units of keV, cm~^s“^, 10“^ ergs~^ cm“^ and lO^'^ergs”^, respectively. 

^ Errors are 68% confidence errors. 


Table 4. Results from the Swift + RXTE spectral analysis^. 



Obs. 5 

Obs. 6 

Obs. 7 

Obs. 8 

Obs. 9 

Obs. 10 

NH(xl022)(cm-2) 

1.01±0.08 

0.74±0.07 

0.84±0.07 

0.80±0.13 

0.77±0.07 

0.93±0.12 

E(l)g.bs (keV) 

4.1±1.6 

7.8±0.7 

8.5±0.5 

7.0T1.0 

7.2±0.6 

8.86±0.20 

cr(l) (keV) 


4.3±0.7 

3.0T0.8 

5.1±1.4 

5.5±0.6 

7±3 

Strength(l) 

181“ 

12±4 

-*-•^-0.4 

28±14 

15±5 

60±30 

E(2)g,bs (keV) 

30.8±0.8 

29.3±1.1 

26.8±0.9 

30.0±1.4 

30.0±2.0 

27.5±2.3 

o-(2) (keV) 

4.7T0.5 

3.4T0.4 

2.7±0.5 

4.0T0.7 

3.5T0.5 

2.15±1.0 

Strength(2) 

19±4 

14±4 

6.0T2.0 

12±4 

16±9 

10±8 

Egauss (keV) 

6.56T0.16 

6.23±0.20 

6.23±0.19 

6.15±0.4 

6.5±0.5 

6.5±0.8 

a (keV) 

U.UU^_Q QQ3 

0.4T0.3 

0.34T0.20 

1.25±0.25 

0.003±0.003 

1.4±0.9 

N(cm-2s-i) 

0.55±0.10 

0.007±0.004 

0.0013±0.0010 

0.05±0.04 

0.0012±0.0007 

0.07±0.05 

r 

A pow. 

-0.02±0.21 

-0.27±0.12 

0.05±0.10 

-0.20±0.20 

-0.49±0.10 

0.0±0.4 

Ecutoff (keV) 

5.4T0.4 

5.8T0.7 

10.3±1.6 

6.0T0.5 

5.2±0.7 

7±4 

Fx ^ (xl0~®ergs“^cm“^) 

4.3±0.4 

6.5±0.6 

6.3±0.6 

6.2±0.6 

3.3±0.3 

3.2±0.3 

Lx(x 103'^ ergs-i) 

2.5±0.3 

3.8±0.4 

3.8±0.4 

3.7±0.4 

2.0±0.2 

1.9±0.3 

CpCA 

1 

1 

1 

1 

1 

1 

Chexte 

1.43T0.20 

0.45±0.15 

0.11T0.05 

0.33±0.06 

0.61T0.15 

0.26±0.17 

CxRT 

0.87T0.08 

0.39±0.03 

0.39±0.03 

0.29T0.03 

0.41T0.03 

0.48±0.05 


0.93(182/195) 

0.97(152/156) 

1.13(183/161) 

0.87(97/112) 

0.94 (97/103) 

1.02 (99/97) 


^ Model used: phabsx constant XgabsXgabs ( cutoffpl + gaussian ) for the 1-60 keV Swift + RXTE spectra. 
^ Errors are 68% confidence errors. 


The fits and the residuals of the whole spectra improved 
(y^/i/Ril.5 — 6) by the addition of a cyclotron line scat¬ 
tering feature (CRSF; hereafter referred to as “cyclotron 
feature”). It was accounted for with a gaussian absorp¬ 
tion line (gabs). The use of this model component and/or 


studies (iKrevkenbohm et al.l. 2005t 

Tsvga 

nkov et al.l. 12006 

R.eie et al.l. 20061: Reis & Nesoolil. 

20131: 

Lutovinov et al. 

201511. The parameters were let to be free in the HEXTE 


data and tied to those in the XRT and PCA data. 


The fits still show significant residuals at 7-10 keV. This 
feature is known as the “10 ke V feature”, and its origin is 
uncertain (ICobiirn et al.l . 120?)^ . We modeled it as a gaus¬ 
sian absorption line (gabs), thus improving the fits substan¬ 
tially (x^/i^wl.O — 1.2). In this case the parameters were 
let to be free in the PCA data and tied to those in the 
XRT and HEX TE data. The value s obtained are within 
those expected ijCoburn et al.l . l2002tl but with broader dis- 
persion than those obtained e.g. for XTE J1946-I-274 by 
iMiiller et al.l (l2012ll . 
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Table 5. Results from the Suzaku spectral analysis^. 


Obs. 14 

Nh(x1022)(j.jjj-2) 

1.2+0.1 

kTbb (keV) 

2.4+0.9 

Nbb 

0.011+0.007 

r 

pow. 

-1.3+1.6 

Ecutoff (k6V) 

15 - 190 

Ef (keV) 

l.O+l.O 

Fx (xl0“^^) (ergs“^ cm“^) 

2.4+0.9 

Lx(xl03®)(ergs-i) 

1.4+0.5 

Cxisi 

1 

CxiS3 

1.35+0.04 

CpiN 

0.4+0.3 


1.35 (221/163) 


^ Model used: phabsx constant xhighecutCdiskbb + 
powerlaw) for the 0.8-30 keV Suzaku spectrum. 

^ Errors are 68% confidence errors. 


5 10 60 1 10 60 

Energy (keV) Energy (keV) 

Fig. 6. The two (upper) panels (left and right) show two 
example spectra of V 0332+53. Obs. 2 is one of the two 
observations with the highest count rates while Obs. 6 is 
the dataset with the highest count rate with all the instru¬ 
ments available (black, red and green correspond to PCA, 
HEXTE and XRT, respectively). The (lower) panels show 
the behaviour of the residuals when adding the different 
components one by one. Number at the right of each panel 
indicate the best fit statistics (x^/d.o.f.). From upper to 
lower panels: 1) best-fit using the cutoffpl model compo¬ 
nent only (i.e. continuum model); 2) residuals after adding 
the RilOkeV (absorption) feature to the continuum model; 
3) residuals after adding the cyclotron feature; 4) residuals 
after adding the Fe line to the model (and final residu¬ 
als of the model). Note that the same vertical axis is used 
in both panels. 


Although the beneath of the spectral fits, there are still 
peaked positive residuals in the region 6-7 keV of the PCA 
spectra. We accounted for them adding a gaussian line pro¬ 
file (gaussian) at 6.4-6.97keV in order to account for the 
Fe Kq fluorescence line. The PCA residuals flattened and 
we obtained the hnal fitting solution in this way. 

To account for uncertainties in the absolute flux normal¬ 
ization between PCA and HEXTE we introduced a multi¬ 
plicative constant that was fixed to 1 for the PCA and let 
to vary freely for HEXTE. In the case of XRT, PCA and 
HEXTE data we fixed the multiplicative constant of PCA 
to 1 and the rest to vary freely. The most relevant results 
of this spectral analysis and the derived unabsorbed fluxes 
and luminosities are shown in Tab. I3l4l and in Fig. [51 

In addition, we fitted the spectra obtained with Suzaku, 
during much a lower luminosity period (see Fig. [7] and 
Tab.[3). Spectra from XISl, XIS3 and PIN were fitted in the 
0.8-10,14-25keV energy range, with the 1.6-2.3keV energy 
range ignored due to the uncertainties in calibration asso¬ 
ciated with the instrumental Si K edg4IE|. The spectra were 
grouped to have at least 50,250 counts for the two XIS and 


PIN, respectively. In this case a model composed by pho¬ 
toelectric absorption (phabs) and a power-law with an ex¬ 
ponential high-energy cutoff of the type “highecut” was a 
rather better description than “cutoffpl” (i.e. = 15 

for the same number of d.o.f and better residuals). With 
this model we obtained a statistics of = 1.4 and 

systematic residuals at E < 4 KeV, that we improved by 
adding a blackbody spectrum with normalization propor¬ 
tional to the surface area (bbodyrad), obtaining a fit statis¬ 
tics of {x^liz = 1.3 for the same number of d.o.f.). Because 
the normalization of the black-body component was not 
well constrained we tried to improve it and changed it by an 
emission component describing an accretion disc consisting 
of multiple blackbody compo nents Idiskbbi lMitsuda et al.l 
11984 iMakishima et al.lIl986ll . We obtained a very similar 
fit statistics (see Tab. 0. A plausible origin for the first 
model component (bbodyrad) is the NS surface or perhaps 
the boundary layer (if present). Low-level accretion onto 
the surface of a neu tron star can indeed produce a black¬ 
body like spectrum (IZampieri et al.l . Il995ll . Still, a disc ori¬ 
gin cannot be ruled out since it was also possible to fit 
the data with a multicolor disc blackbody. Similar results 
for a set of low-accre t ing N S have been recently obtained 
(jArmas Padilla et al.l . l2013ll . but the temperatures found 
are kT« = 0.5 — 0.7keV, thus much lower than in our case 
(kTsi = 2.4±0.9keV). The higher value we obtain might 
be indicative of the much higher accretion rate in our case. 
Nevertheless, with the current data it is not possible to 
discern which model component is better to describe the 
data. 

To account for the uncertainties in the absolute flux 
normalization between XISl, XIS3 and PIN we introduced 
a multiplicative constant that was fixed to 1 for XISl and 
let to vary freely for XIS3 and PIN. The most relevant 
results of this spectral analysis and the derived unabsorbed 
fluxes and luminosities are shown in Tab. [5] 

4. Discussion 

4.1. Long-term Be-disc/neutron star interaction 


From the long-term observations of V 0332+53 we can see 
http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/sical.html that V 0332+53 spends most of its time in the optical 
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OBS. 14 (SUZAKU) 



ENERGY (keV) 

Fig. 7. Energy spectrum from V 0332+53 during Obs. 14 
with Suzaku (see Tab. [S] for the observation details and 
the results of the spectral analysis). Black, red and green 
correspond to XISl, XIS3 and PIN, respectively. 



Fig. 8. Flux in the 5-60 keV energy range versus the cy¬ 
clotron line energy centre for Obs. 1-11. See on-line material 
Tab. [3] for the (not-rounded) list of values used. 


brightening phases (Fig [3]). The most significant one is the 
period between 1993-2005. It started approx. 4 years after 
the giant outburst that the system exhibited in 1989 and 
peaked simultaneously with the X-ray maximum in 2004. 
Indeed all of the optical brightening episodes of BQ Cam 
are accompanied with the X-ray activities. It is important 
to note that the brightness of the optical companion de¬ 
creases after the X-ray outbursts. This behavior can be well 
explained by the weakening of the decretion disc during the 
mass transfer to the NS. 

Having a moderate eccentricity of e = 0.31 with a short 
orbital period of Porb = 34.25 d, V 0332+53 is one of the 
BeXBs for which the effect of the truncation exerted by 
the NS is likely to be observed. Due to the small orbit of 
the companion, the decretion disc cannot expand so much. 
Thus, it is truncated at a radius smaller than the critical 
Roche radius. Since the amount of the material supported 
by the Be star to the disc would be accumulated in time, the 
increase in the optical brightness can be understood in this 
way. However the accumulation of the material particularly 
at the outer part of the disc makes these regions to become 
optically thick. According to theory by either the radiation¬ 
driving warping or a global den sity wave the outer part 
of th e disc is strongly elongated ([Okazaki and Negueruelal . 
l200lh . In general, the existence of the global density waves 
in the decretion disc are observationally supported by the 
V/R variations seen in the emission line profiles. Although 
the line profile variations of V 0 332+53 were previously 
observed ([Neeueruela et al.l . 1199811 during the period 1990- 
1991, Hq, lines are always seen in nearly symmetric single- 
peaked emissions without any significant variations for the 
time interval 2006-2014. Therefore, we do not have any ev¬ 
idence to support the idea of the perturbations occurred in 
the disc as well as the V/R variations. 

4.2. X-ray behaviour 
4.2.1. Spectral evolution 

From the RXTE and Swift observations reported in 
this work we see a decrease of the X-ray luminosity of 


V 0332+53 during the latest part of the intermediate- 
luminosity 2008 outbu rst. Assuming a 7kpc distance 
(iNegueruela et al.Lll999ll . the luminosity gradually dropped 
from ^6 to ^IxlO^^ergs”^ in the 5-60 keV energy range 
from Obs. 1-11. The luminosities are also in accordance 
with t hose obtained during the decay phase of the 2005 out¬ 
burst (iMowlavi et al.l . l2006ll . Eventually, during the Suzaku 
observation (Obs. 14) we obtain the lowest luminosity 
(~lxl0^® ergs“^). 

A recent study of V 0332+53 during giant outbursts 
at different luminosities has provided cru cial insights into 
the m ost frequent states of V 0332+53 (iReig fc Nesooli 
I 2 OI 3 II . Their Fig. 2 shows the existence of two separate 
regimes (the high-luminosity diagonal branch and the low- 
luminosity horizontal branch). They are separated by a 
critical luminosity (Lcrit~l — 4xl0^^ergs“^), which cor¬ 
responds to the range of luminosities reported in our 
work during Obs. 1-11. Both regimes/states are clearly 
well-populated by X-ray observations whilst observations 
around the critical luminosity are scarce. 

The spectral variability of V 0332+53 has been well 
studied by different authors, confir ming the softe n ing o f 

V 0332+53 as the flux decreased. iBavkal et al.l (1200711 
suggested that this type of spectral softening with de¬ 
creasing flux was mainly a c onsequence of m ass accre¬ 
tion rate change. In particular, iReig et id] (|20inll (and ref¬ 
erences therein) explained that increased mass accretion 
rates are expected to result in harder X-ray spectra due 
to Comptonization processes. In our work we see the same 
trend photon index versus flux. 

The values previously obtained for the energy centre 
of the cyclotron line have shown a stro ng anti-correlation 
with luminosity in V _ 033 2+53 (iTsvgankov et al 


20061: iMowlavi et al l . 120061 : ITsvgankov et al.l . I201C : 


Reig fc Nesooli . 201311 during giant outbursts. Somehow 


weaker anti-correlat i on has been found in 4U 0115+63 
(iMihara et al.l . l2004 iNakaiim a et ah. 200611 and no cor- 
relation in lA 0535+262 ([Caballero et ah . 200711 . In our 
work a positive correlation of the line energy centre 
with luminosity is observed (see Fig. [S]). The correlation 
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coefficient of the fit is +0.64 (i.e. null two-tailed probability 
of 0.03). This kind of behaviour has also been observed in 
A 0535+262 during an intermediate-lum inosity outburst 
(Lx = 4x l0^^ergs~^:[^ rtore et al.ll2015ll . We notice that 
previously iReig fc Nespolil ( 201311 already pointed out for 
a tentative positive correlation for the latter source. This 
effect appears to be a natural consequence of accretion 
onto X-ray pulsars with the strongest magnetic fields 
in the sub-critical regime (|Mushtukov et all . l2015f) . The 
values we have obtained for the width of the cyclotron line 
are all w4keV, in contrast to the values obtained during 
high flux states (6-8 keV; iTsveankov et al.l l2010ll but in 
full agreemen t with those obtain ed in intermediate flux 
states 14 keV: iMowlavi et al.ll2006ll . 


4.2.2. Quasi Periodic Oscillations 


We have found that the PDS of V 0332+53 during our 
observations are best described by band-limited noise plus 
QPO noise during the latest observations. The QPO noise 
consists of a component at a frequency of ssO.22 — 0.23 Hz. 
This feature has also been reported in the 2-60 keV RXTE- 
PCA PDSs obta ined during the 2004-2005 outburst decay 
(|Qu et al.l . l2005fl . They found a QPO at «0.22Hz in addi¬ 
tion to a narrow peak (i.e. pulsation from the pulsar) close 
to that frequency (Ri0.23Hz). The QPO was detected when 
the luminosity was at a mid-flux or even lower flux level of 
the outburst. The rms strength of this QPO also amounts 
to a few per cent in that work. The 0.22 Hz QPO was ob- 
served around the sp in frequency in both the INTEGRAL 
(|Mowlavi et ah . ^00611 and RXTEdata dur i ng the 2005 out- 
burs t decline (iQn et al.l. 120051 : iReig et rill . l2006ll . In addi¬ 
tion, [^keshi^^^^OCO^^discovered another QPO with 
a centroid frequency of RS0.05HZ and a relative root mean 
square (rms) amplitude of ~5%, in the 2.3-37.2 keV data 
from the 1989 outburst of V 0332+53. We do not find any 
signal of the QPO at ~0.05 Hz although we find the pres¬ 
ence of a broad component at around that frequency (that 
may turn eventually into peak noise, i.e. QPO). The QPO 
detected in V 0332+53 at ss0.22Hz constitutes the first de¬ 
tection of a QPO riding on the spin frequency of a NS and 
points out to a strong coupling between the periodic (i.e. 
the pulses) a nd the red-noise (i.e. broad band QPO signal) 
components (|On et al.l . l2005ll . 

The existence of a QPO signal with a similar frequency 
to the spin of the NS suggests that they are related some- 
ho w. One model t o expl ain s uch coinciden c e was given 
bv iLazzati fc Stellal (Il997ll and iBnrderi et al.l (Il997ll . They 
proposed that the QPO is originated in the magnetosphere 
around the NS. Additionally, they argued that the accre¬ 
tion flow is inhomogeneous near the surface of the NS, and 
that random shots are characterized by an arbitrary degree 
of modulation. The combination of rotation and radiative 
transfer effects should therefore produce a periodic modu¬ 
lation of the shots similar to that of any continuum X-ray 
emission from the polar caps. They proposed that a cou¬ 
pling between the periodic and red-noise variability should 
be frequently present in X-ray pulsars. 

In our case an d in contras t to pr e vious observ ations dur¬ 
ing giant flares (IReig et al.l . 120061 : iReid . 1200811 the QPO 
signal appears during the lowest of the measured fluxes 
(Lx<4xI0^^ ergs“^). If the QPO has an origin on the mag¬ 
netosphere of the NS (and is not an artifact from our ob¬ 
servations) , it should be “more visible” when the accretion 


column is more tenuous (or even absent). Summarizing, it 
looks like that during our latest Obs. 6-II we might be wit¬ 
nessing the innermost parts closest to the magnetosphere 
surroundings. 


5. Summary and conclusions 

In this paper we present a multi-wavelength study of the 
Be/X-ray binary system V 0332+53 with the main goal 
of studying its transient behavior during its intermediate 
luminosity X-ray outburst in 2008. After being R::!3.5yr in 
quiescence in X-rays the system underwent a new active pe¬ 
riod. A new X-ray outburst was detected around 17 October 
2008, roughly six days before the periastron passage of the 
NS. We note that the optical companion entered the bright¬ 
ening state nearly 1.5 yr before the X-ray outburst. 

V 0332+53 was in quiescence until the optical magni¬ 
tude approached a peak value of V = 15.2 mag. at the end 
of 2009, indicating the beginning of another mass ejection 
event. In November 2009 a series of five small Type I peri¬ 
odic outbursts was observed. The X-ray activities finished 
by the end of May 2010 whist the optical magnitude was 
still fading. It ceased around August 2011 for about Rs3.5yr. 
Our recent data showed that a mass ejection event is cur¬ 
rently taking place, exhibiting an increase of ssO.6 in the V 
mag. It probably started around 2013 and peaked in June 
2015 when a new (type H) X-ray outburst was detected. 
This is similar to those observed during the brightening 
episodes occurred e.g. in 1983 and in 2004, but higher than 
the ones in 2008 and 2009. Our only IR measurement from 
December 2014 also showed an unusually bright Be com¬ 
panion. 

The broad-band 1-60 keV X-ray spectra of the NS dur¬ 
ing the decay of the 2008 outburst were all well fitted 
with the standard cutoffpl phenomenological model, en¬ 
hanced by a narrow iron fluorescence line at 6.4 keV to¬ 
gether with a cyclotron line scattering feature at ss30 keV. 
Our study confirms the softening of V 0332+53 as the 
flux decreased. During our lowest flux observation made 
with Suzaku (Lx = 10^®ergs“^) on 2010, compatible with 
V 0332+53 being at a level slightly above from quiescence, 
we detected a very soft spectrum that could be well de¬ 
scribed by adding a (single or multiple) soft black-body 
component (s). 

We tentatively see an increase of the cyclotron line en¬ 
ergy with increasing flux. If confirmed (with better obser¬ 
vations) this might constitute the first detection of such 
positive correlation of the cyclotron line energy versus lu¬ 
minosity in V 0332+53. Regarding the fast aperiodic tim¬ 
ing properties, we detect a QPO at 227±9mHz during the 
lowest luminosities. The latter might indicate that the in¬ 
ner regions surrounding the magnetosphere are more visible 
during the lowest flux stages. Due to their low significance 
we advice that longer and better observations with current 
(and/or future) X-ray satellites are needed in order to con¬ 
firm/discard these results. 
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Table 1. ON-LINE MATERIAL. Optical photometric 
observations from the IAC80 telescope. 


M,JD 

V 

error 

B 

error 

56871.2368550198 

15.460 

0.013 

17.037 

0.071 

56993.9640124198 

15.360 

0.006 

16.919 

0.021 

57009.0723979599 

15.373 

0.005 

17.006 

0.013 

57026.9773657001 

15.372 

0.008 

16.945 

0.031 

57036.8705536202 

15.373 

0.005 

16.967 

0.013 

57054.0418117600 

15.332 

0.008 

16.853 

0.031 

57054.8470330602 

15.364 

0.006 

16.955 

0.024 
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Table 2. ON-LINE MATERIAL. Hq line EW and FWHM measurements of V 0332+53. Negative values of EW 
indicate that the line is in emission. 


DATE 

MJD 

EW 

(A) 

FWHM 

(A) 

Telescope 

26-Sep-2006 

54004.9238 

-6.21+0.71 

6.18+0.38 

RTT150 

23-Dec-2006 

54092.9543 

-7.05+0.56 

5.83+0.55 

RTT150 

07-Apr-2007 

54197.7388 

-5.85+0.33 

6.48+0.83 

RTT150 

19-Jul-2007 

54300.0279 

-5.11+0.10 

6.83+0.22 

RTT150 

20-Jul-2007 

54301.0591 

-4.80+0.49 

6.40+0.68 

RTT150 

13-Sep-2007 

54356.9962 

-3.82+0.18 

6.20+0.67 

RTT150 

06-0ct-2007 

54379.0510 

-5.26+0.18 

6.62+0.40 

RTT150 

14-Dec-2007 

54448.8768 

-4.62+0.17 

6.42+0.29 

RTT150 

26-Aug-2008 

54704.9757 

-4.70+0.09 

6.31+0.48 

RTT150 

22-Sep-2008 

54731.0652 

-3.97+0.09 

6.49+0.62 

RTT150 

29-Oct-2008 

54768.8500 

-4.45+0.10 

6.69+0.26 

RTT150 

04-Mar-2009 

54894.7260 

-3.13+0.23 

6.36+0.34 

RTT150 

24-Sep-2009 

55098.9814 

-5.21+0.12 

6.56+0.44 

RTT150 

21-Oct-2011 

55855.9805 

-6.44+0.10 

7.84+0.19 

OSN-150 

16-Jan-2012 

55942.8952 

-4.38+0.07 

7.22+0.28 

RTT150 

17-Jan-2012 

55943.7739 

-4.86+0.07 

9.38+0.62 

RTT150 

31-Jan-2012 

55957.8646 

-4.08+0.11 

7.97+0.48 

OSN-150 

31-Jan-2012 

55957.8864 

-4.78+0.22 

9.15+0.49 

OSN-150 

05-Feb-2012 

55962.8764 

-5.45+0.10 

7.42+0.20 

RTT150 

04-Mar-2012 

55990.8104 

-4.22+0.09 

6.82+0.17 

RTT150 

27-Mar-2012 

56013.8502 

-3.61+0.23 

7.78+0.45 

OSN-150 

27-Mar-2012 

56013.8618 

-4.31+0.31 

8.49+0.49 

OSN-150 

29-Jul-2012 

56137.2382 

-4.84+0.30 

7.33+0.57 

NOT 

08-Sep-2012 

56178.1507 

-6.56+0.32 

10.05+0.06 

OSN-150 

09-Sep-2012 

56179.2028 

-4.79+0.20 

6.77+0.57 

NOT 

16-Sep-2012 

56186.0446 

-5.83+0.10 

6.84+0.28 

RTT150 

lO-Oct-2012 

56210.1002 

-6.14+0.18 

7.01+0.03 

OSN-150 

09-NOV-2013 

56605.0347 

-7.44+0.08 

13.73+0.66 

RTT150 

05-Jan-2014 

56662.8496 

-5.25+0.14 

6.48+0.67 

RTT150 

05-Feb-2014 

56693.7630 

-4.34+0.20 

6.62+0.37 

RTT150 

06-Feb-2014 

56694.8018 

-4.93+0.33 

7.28+0.40 

RTT150 

20-0ct-2014 

56951.0575 

-4.92+0.07 

6.96+0.48 

RTT150 

21-Oct-2014 

56951.0563 

-4.81+0.11 

6.32+0.13 

RTT150 

22-NOV-2014 

56983.9104 

-5.23+0.13 

7.27+0.20 

RTT150 
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Table 3. ON-LINE MATERIAL. Not-rounded values 
from Tab. [3] used in Fig. |8] 


Obs. 

E(2),,bs (keV) Fx(ergs-icm-2) 

1 

31.6+0.7 

(1.0+0.4)xl0"® 

2 

31.3+0.8 

(0.8+0.4)xl0"® 

3 

31.8+0.8 

(4.70+0.10) X 10”'^ 

4 

29.7+1.2 

(4.5+0.3)xl0"® 

5 

30.8+0.8 

(2.7+1.5)xl0"® 

6 

29.3+1.1 

(3.2+2.4)xl0"® 

7 

26.8+0.9 

(3.1+2.5)xl0"® 

8 

30.0+1.4 

(3.0+2.5)xl0"® 

9 

30.0+2.0 

(1.7+1.3)xl0"® 

10 

27.5+2.3 

(1.5+1.3)xl0-® 

11 

28.4+1.0 

(1.3+0.7)xl0"® 


^ Fluxes are in the 5-60 keV energy range. 
^ Errors are 68% confidence errors. 
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